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ABSTRACT 

We analyze observations from the Interface Region Imaging Spectrograph (IRIS) of the Mg Ilk line, 
the Mg II UV subordinate lines, and the O I 135.6 nm line to better understand the solar plage chro¬ 
mosphere. We also make comparisons with observations from the Swedish 1-m Solar Telescope (SST) 
of the Ha line, the Ca II 8542 line and Solar Dynamics Observatory/Atmospheric Imaging Assembly 
observations of the coronal 19.3 nm line. To understand the observed Mg II profiles, we compare these 
observations to the results of numerical experiments. The single-peaked or flat-topped Mg II k profiles 
found in plage imply a transition region at a high column mass and a hot and dense chromosphere 
of about 6500 K. This scenario is supported by the observed large-scale correlation between moss 
brightness and filled-in profiles with very little or absent self-reversal. The large wing width found in 
plage also implies a hot and dense chromosphere with a steep chromospheric temperature rise. The 
absence of emission in the Mg II subordinate lines constrain the chromospheric temperature and the 
height of the temperature rise while the width of the O I 135.6 nm line sets a limit to the non-thermal 
velocities to around 7 km s“^. 

Subject headings: Sun: atmosphere — Sun: chromosphere — Sun: faculae, plages 


1. INTRODUCTION 

Plages are regions of the solar atmosphere with strong 
unipolar magnetic field. The Interface Region Imaging 
Spectrograph mission (IRIS; De Pontieu et al. 2014b), 
provides new diagnostics for the properties of plage in 
the form of the Mg II h&k lines, the Mg II subordi¬ 
nate UV lines, which are sensitive to heating in the low 
chromosphere (Pereira et al. 2015); and the O I line at 
135.56 nm that samples the chromosphere but forms un¬ 
der optically thin conditions (Lin & Carlsson 2015). In 
this Letter we investigate what the IRIS observations, 
together with observations from the Swedish 1-m Solar 
Telescope (SST; Scharmer et al. 2003) and the Atmo¬ 
spheric Imaging Assembly (AIA; Lemen et al. 2012) on 
board the Solar Dynamics Observatory (SDO) can tell 
us about the heating of the chromosphere and the neigh¬ 
boring regions (transition region (TR) and corona). 

Understanding the physical mechanisms driving the 
formation of Mg II h&k is of great interest, not only 
for understanding the heating in the chromosphere, but 
also because these lines are used extensively as a proxy 
for the Sun’s cyclical variability over the entire UV spec¬ 
tral domain (i.e., the “Mg II index”; Heath & Schlesinger 
1986). 

2. OBSERVATIONS 

We used four data sets of active regions (including 
sunspots, pores, and plage regions) obtained with IRIS, 
one of which (2014 June 11 at 07:36 UT) was coordi¬ 
nated with Ha and Ca II 854.2 nm spectral line scans 
from the CRisp Imaging SPectrapolarimeter (CRISP; 
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Scharmer 2006; Scharmer et al. 2008) mounted on the 
SST. All IRIS raster scans were dense (0?35 steps) with a 
spatial sampling of 0?17 along the slit (unless otherwise 
noted). The 2014 June 11 data set focused on NO A A 
active region (AR) 12080 and consists of 20 large dense 
rasters with 96 raster steps resulting in a field of view 
of 3375x182" centered at {x,y) = (57375,-19973). The 
exposure time per raster step was 4 s. The raster cadence 
was 516 s. 

The other IRIS data sets are similar large dense rasters 
with 400 raster steps and a resulting field of view of 
140" xl82". They have different exposure times, raster 
durations, and pointings of, respectively, 30s, 3.5 hr, 
{x,y) = (-113"2,-23470) (NOAA AR 12104) for the 
2014 July 4, 11:40 UT data set; 8s, 1 hr, (x^y) = 
(-108", 106") (NOAA AR 12139) for the 2014 Au¬ 
gust 16, 21:17 UT data set; and 30s, 3.5 hr, {x,y) = 
(1373, -24374) (NOAA AR 12187) for the 2014 October 
18, 06:45 UT data set. The last data set included spatial 
binning by two along the slit to improve the signal-to- 
noise of the weak O I 135.6 nm line. All data sets were 
co-aligned with images in the 17.1 and 19.3 nm AIA pass- 
bands. 

3. OBSERVED PLAGE PROPERTIES 
3.1. Typical Mg Ilk Profiles 

The Mg II k line shows a wide range of profiles through¬ 
out AR 12104 (Fig. 1). While the average profile of the 
field of view shows the typical central reversal (although 
significantly reduced compared to the quiet Sun), this av¬ 
erage hides a wide range of types of profiles. Quiet-Sun 
regions in the vicinity of a plage (which is bright in Mg II 
wing; Fig. 1(a)) show deep central reversals with typi¬ 
cally fainter k 2 peaks (for a definition of spectral features, 
see Leenaarts et al. 2013a) and mostly narrower profiles 
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Fig. 1. — IRIS observations of AR I2I04 in the wing of the Mg Ilk line (a), at the Mg Ilk core (ka) (b) and in the subordinate Mg II 
blend at 279.88 nm (c). Typical plage Mg Ilk profiles are shown in (d) for the diamond locations in (a); with non-plage Mg Ilk profiles in 
(e) for the star locations in (b). Sample plage profiles of the subordinate Mg II 279.88 nm blend are in (f) for the same diamond locations 
as in (a). Full field of view average profile dashed in panels (d)-(f). Light-blue contours in (c) show where the subordinate blend is in 
emission, the full black line in (f) is an example profile at the black square location in (c). 


(blue, cyan, and green profiles/locations in Fig. l(b),(e)). 
Sunspots (green in Fig. l(a),(d)) and pores (yellow in 
Fig. 1(b), (e)) are fainter, much narrower than the aver¬ 
age AR profile (dashed lines in Fig. l(d)-(f)) and typi¬ 
cally show a single-peak profile. Profiles in plage are very 
different from quiet-Sun and sunspot/pore profiles: they 
are brighter, wider, and typically either single-peak (yel¬ 
low, Fig. 1(a), (d)), fiat-topped (blue. Fig. 1(a), (d)) or 
with a very small central reversal (cyan. Fig. l(a),(d)). 

Figure 1(c) shows that Mg II 279.88 nm is indeed sen¬ 
sitive to the chromosphere with many morphological fea¬ 
tures similar to those seen in Mg II ks (b). While some 
locations show this subordinate blend in emission (light 
blue contours in (c); black square/profile in (c)/(f)), they 
are typically in sunspots, explosive events such as bombs 
(Peter et al. 2014) or locations in the canopy that sur¬ 
round the plage. The subordinate blend almost never 


goes into emission in the plage regions. This appears to 
provide strong constraints on the chromospheric temper¬ 
ature profile (see § 4). 

While plage profiles without a central reversal are 
common, not all plage regions show such profiles. On 
large, active-region-size, spatial scales, most plage re¬ 
gions are bright in k 2 (defined as average brightness in 
two fixed wavelength ranges to the blue/red of the line- 
center wavelength, as derived from the average profile; 
Fig 2(a)). On the same scales, the ks map (brightness 
at the wavelength of line center of the average profile; 
Fig 2(b)) shows more variability with some bright k 2 
regions associated with fainter ks intensity, and others 
showing very bright ks- This is illustrated by panel (c): 
for single-peak profiles ks intensity is larger than k 2 , so 
regions with many single-peak profiles appear dark in 
panel (c). In summary, most plage regions are bright in 
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Fig. 2. — Correlations and comparisons between various Mg II k profile parameters and SST and SDO/AlA observables. Top row: Mg II k 2 
(a) and Mg Ilka (b) of AR 12139, the difference between k 2 and ka (c) and AIA 19.3 nm (d). Bottom row: Mg Ilka intensity of a small 
plage region of AR 12080 (IRIS, e), Ha line core (SST, (f)), Ca II 854.2 line core (SST, (g)) and AIA 19.3 nm (SDO, (h)). This figure is 
accompanied by two online animations. 


k 2 , but only a subset show more single-peak profiles. 

3.2. Correlations with the TR and eorona on large 

seales 

Comparison with the SDO/AIA 19.3 passband (dom¬ 
inated by Fe XII, logT^ax ~ 6.2; Fig. 2(d)) shows that 
such single-peak regions have a spatial correlation, on 
large spatial scales, with regions of enhanced moss emis¬ 
sion. This AIA passband contains both bright coro¬ 
nal loops, dark absorbing features, and bright, low-lying 
moss (e.g., (100", 95")). Moss regions are the upper 
TR footpoints of high pressure coronal loops that are so 
hot that their TR emission is formed at 1 MK, so that 
the SDO/AIA 19.3 nm emission originates from low-lying 
plasma that occurs at the same heights as chromospheric 
jets or dynamic fibrils (Berger et al. 1999; Fletcher & De 
Pontieu 1999; De Pontieu et al. 1999) that absorb some 
of the EUV emission and give the mottled appearance 
typical for moss regions (see, e.g., bottom row of Fig. 2). 

The brightness of moss is thought to be a good proxy 
for and linearly related to the pressure in the overlying 
hot coronal loops (Martens et al. 2000). We thus find 
that the brightest moss regions, i.e., locations with high 


coronal pressure, and thus TR at higher column mass, 
typically show more single-peak profiles. 

3.3. Correlations with ehromosphere and TR on small 

seales 

We also find correlations on the smallest observable 
spatial scales (bottom row of Fig. 2). We find that 
on sub-arcsecond scales the brightness in a small mossy 
plage region of Mg IIks, Ha line center, and Ca II 854.2 
nm are well correlated. This correspondence is in part 
due to the small dark features visible in all three lines 
that are associated with dynamic fibrils (Hansteen et al. 
2006; De Pontieu et al. 2007). This indicates that these 
diagnostics are all sensitive to upper chromospheric con¬ 
ditions, including velocity fields. The Ca II 854.2 nm 
brightness and Ha width have been proposed as a proxy 
for chromospheric temperatures in more quiescent condi¬ 
tions (Cauzzi et al. 2009). This relationship seems doubt¬ 
ful in plage: (1) the Ha width is not well correlated with 
Ca II 854.2 nm brightness (see the animation of Fig. 2); 
(2) Ca II 854.2 nm is well correlated with Mg II ks bright¬ 
ness, which as shown below appears to be more sensitive 
to the TR rather than chromospheric conditions. 
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At small spatial scales, there is also a good, though not 
perfect, correspondence with AIA 19.3 nm moss emis¬ 
sion, especially with dark moss features. This is perhaps 
not surprising as these features have been previously as¬ 
sociated with bound-free absorption from dynamic fibrils 
(De Pontieu et al. 1999). On the other hand, many of 
the bright moss locations have an equivalent bright re¬ 
gion in Mg IIks, Ha line center, and Ca II 854.2 nm line 
center. These findings expand on the previously found 
relationship between upper TR moss emission and Ha 
line center (De Pontieu et al. 2003) or Lyman-a emission 
(Vourlidas et al. 2001). 

3.4. Properties of Mg II k in Plage 

Figure 3 shows various observables of AR 12187. Panel 
(a) shows context by showing the intensity at 280 nm, 
a wavelength between the Mg II h&k lines where the 
intensity is formed in the upper photosphere. Magnetic 
areas are clearly seen as increased intensity, and there is 
a large plage area in the left part. We use this intensity 
to determine a mask for plage (green contours). 

The width of the O I 135.6 nm line (panel (b)) in the 
plage area is remarkably constant, around 7.8 km s“^ 
with a small spread (panel (g)), and is larger than in 
internetwork areas. The O I 135.6 nm line is optically 
thin, and the width gives a direct measure of the width 
of the atomic absorption profile (Lin & Carlsson 2015). 
At a temperature of 7kK the thermal 1/e width of the 
O I 135.6 nm line is 2.7km s“^ leading to a mean non- 
thermal width of 7.3 km s“^. 

The wing width of the Mg Ilk line (panel (c)) is simi¬ 
larly larger in the plage area than in non-plage areas with 
a mean 1/e width of 30 km s“^ with a small spread (g). 
The large wing width of the Mg Ilk line comes from a 
large ’’opacity broadening factor” (see Rathore & Carls¬ 
son (2015) for a discussion). The wings outside the k 2 
peaks sample the velocity field in a similar region as the 
O I 135.6 nm line. There is a clear correlation between 
the widths of the two lines but not one-to-one, showing 
that the Mg II k wing width is not only influenced by the 
velocity field in the chromosphere but also by a varying 
opacity broadening factor. 

The radiation temperature of the Mg Ilk peak inten¬ 
sity is around 6400 K in plage with a smaller spread 
compared with the full field of view (panels (d),(h)). The 
plage area shows a mottled appearance with very little 
correlation with the intensity at 280 nm on small scales. 
There is some correlation with the Mg Ilk wing width 
(c). 

Panel (e) shows how the central reversal of the 
Mg Ilk line is filled in by showing the difference in ra¬ 
diation temperature between the peak intensity and the 
ks intensity. Single-peak profiles thus have a difference 
of zero. Most of the central part of the plage has a very 
weak central reversal or single-peak profiles. These loca¬ 
tions also have larger maximum intensity (d) and larger 
Mg Ilk wing width (c). Stronger self-reversals occur in 
small patches outside the central part of the plage. These 
patches are correlated with asymmetric peaks, larger 
peak intensity, and larger k 2 peak separation (panel (f)) 
but smaller Mg Ilk wing width (c). 

The Mg II k 2 peak separation (panel (f)) shows a lot 
more structure than the Mg Ilk wing width (c). The 


separation is large for the asymmetric profiles that have 
a rather deep central reversal (bright patches in the non¬ 
central part of the plage in panel (e)). 

4. MODELING 

The Mg II h&k line profiles from plage regions show 
some common features: their emission cores are wider 
than in other atmospheric regions, they are brighter, and 
their central depressions are shallower or absent. 

By trial and error we constructed a model that matches 
the observed average Mg Ilk profile reasonably well by 
simultaneously solving the non-LTE problem for hydro¬ 
gen, calcium, and magnesium, including charge conserva¬ 
tion and enforcing hydrostatic equilibrium using the RH 
code (Uitenbroek 2001). For hydrogen and calcium, we 
used the standard 5-level-plus-continuum H I and Ca II 
models that come with RH; for magnesium, we used the 
10-level-plus-continuum Mg II atom from Leenaarts et al. 
(2013a). The Lya, Ly^, Ca II H&K and Mg II h&k lines 
were computed including partial redistribution; all other 
lines were computed assuming complete redistribution. 

The model is constructed from the photosphere and 
TR of the one-dimensional static FALP model (Fontenla 
et al. 1991). Compared to FALP, our model has an ex¬ 
tended temperature minimum, a steeper chromospheric 
temperature rise, a constant chromospheric temperature 
plateau, and the TR located at a larger column mass. 
In addition, we replaced the FALP non-constant micro¬ 
turbulence in the temperature plateau with a constant 
value. In Figure 4 and 5 we show the main parameters 
of our model and compare the synthetic Mg II k linecore 
and the subordinate blend at 279.88 nm with the average 
observed plage spectrum. 

We do not propose this model as a realistic model of the 
atmospheric structure of plage; it is a numerical experi¬ 
ment exploring the constraints that the Mg II profiles set 
on the structure of plage chromospheres. We therefore 
explored the sensitivity of the Mg II lines to variations of 
our model by varying the column mass of the TR (pan¬ 
els (a)-(c) in Figure 4), the temperature of the chromo¬ 
spheric plateau (panels (d)-(f) in Figure 4), the micro¬ 
turbulence in the chromospheric plateau (panels (a)-(c) 
in Figure 5), and the column mass of the chromospheric 
temperature rise (panels (d)-(f) in Figure 5). 

Varying the location in column mass of the TR has a 
large effect on k 2 and ks: the emission peaks get higher, 
and the ks minimum fills up as the TR is moved toward 
larger column mass, but even the best-fit model still has 
some central reversal. The subordinate blend gets shal¬ 
lower. 

The temperature of the chromospheric plateau has a 
significant effect on the Mg Ilk core-width: larger tem¬ 
peratures mean larger widths. For temperatures above 
6.5 kK the k 2 intensity increases. The subordinate blend 
reacts strongly to changes in temperature. For temper¬ 
atures below 6.5 kK, it is in absorption; for higher tem¬ 
peratures, it gets marked emission peaks. 

Varying the microturbulence in the chromospheric 
temperature plateau changes the width of the 
Mg Ilk emission core, and the width of the absorption 
line caused by the UV subordinate transitions. 

Changing the location of the chromospheric tempera¬ 
ture rise has two effects. First, the deeper in the atmo¬ 
sphere it is located, the wider is the Mg II k emission core. 
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Fig. 3.— AR 12187 observables: (a) Intensity at 280 nm (photospheric context); (b) 1/e width of a single-Gaussian fit to the O I 135.6 nm 
line, (c) 1/e width of single-Gaussian fit to the outer wings of Mg Ilk (core between k 2 peaks excluded), (d) mean intensity of Mg II k 2 
peaks (or maximum intensity for single peaks), (e) difference between peak intensity of the Mg Ilk line and ka central depression, (f) peak 
separation between Mg II k 2 peaks, (g) histogram of 1/e widths of O I 135.6nm (black) and Mg Ilk line wing fit (red) for the full field of 
view (dashed) and plage (solid), (h) histogram of the radiation temperature of the Mg II k 2 peaks for full field of view (dashed) and plage 
area (solid). All intensities are in radiation temperature. Green contours in panels (a)-(f) show the plage region, with ranges for the color 
table in the lower right corner. 
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Fig. 4.— Sensitivity of the emergent Mg Ilk line to variations in a ID static plage atmosphere model. Top row: (a) best-fit model 
atmosphere (red) as function of column mass. Curves in shades of gray show variations of the column mass of the TR. (b) Emergent 
Mg Ilk line-core profile from the best-fit atmosphere (red) and the variations shown in the panel on the left (corresponding shades of gray). 
The average plage profile from the IRIS observations is shown in blue. The outer minima are photospheric lines from different elements 
than Mg that are not included in the modeling, (c) Same as (b), but now for the two subordinate lines at 279.875 nm and 279.882 nm. 
Bottom row: same as the top row, but now for variations in the temperature of the chromospheric plateau. 
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Fig. 5.— Sensitivity of the emergent lines to variations in a ID static plage atmosphere model. This figure follows the same format as 
Figure 4. Top row: variation in chromospheric microturbulence. Bottom row: variation in column mass of the chromospheric temperature 
increase. 
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Second, deeply located temperature rises lead to emission 
in the wings of the subordinate blend. 


5. DISCUSSION AND CONCLUSIONS 

Our results on the peculiar Mg Ilk profiles found in 
plage provide an intriguing picture of heating in the 
plage chromosphere and hotter layers above. Our nu¬ 
merical experiments suggest that single-peaked or flat- 
topped profiles found in plage arise naturally if the chro¬ 
mosphere is hot and dense, thus ensuring coupling of 
the source function to the temperature, and if the col¬ 
umn mass of the TR is high so that k 2 and ks form at 
the base of the TR. The latter occurs naturally at the 
footpoints of hot, dense coronal loops in which a strong 
thermal conductive flux pushes the TR to high column 
mass. This scenario is supported by the observed large- 
scale correlation between moss brightness (a proxy for 
coronal pressure) and filled-in profiles with very little or 
absent self-reversal. The correlation is not perfect since it 
depends on both chromospheric and coronal conditions. 
For example, cooler (^1 — 2 MK) but dense coronal 
loops can also have high coronal pressures (and thus TR 
at high column mass), yet they would lack signatures of 
AIA 19.3 nm “moss”. In addition, locations with TR at 
high column mass may occur even for a chromosphere 
that is lower density or cooler (leading to centrally re¬ 
versed profiles). Despite these caveats, we often do find 
a reasonable correlation on small, arcsecond spatial scales 
between ks brightness and upper TR moss emission. At 
these small scales, the Ha width does not correlate with 
ks brightness, but the Ca II 854.2 nm brightness does. 
This is unlike the quiet-Sun. 

We find that, contrary to quiet Sun profiles (Leenaarts 
et al. 2013a, b), the ks and k 2 properties in plage 
are sensitive to conditions higher up, at the very 
top of the chromosphere. The Mg Ilk sensitivity to 
(mid-)chromospheric conditions can be found in the 
width of the Mg Ilk wing, the (lack of) emission of the 
subordinate blend at 279.88 nm and the k 2 peak sepa¬ 
ration. Compared to quiet Sun, the Mg II k wing width 
is significantly higher (but remarkably constant around 
^ 30 km s“^) in plage regions. Comparison with non- 
thermal line broadening of the O I 135.6 nm line shows 
that it is smaller than 10 km s“^ and that there is a gen¬ 
eral correspondence on large spatial scales, but not on 
small spatial scales. This suggests that both increased 
microturbulence as well as so-called opacity broadening 
must play a role in the large wing widths of Mg II k. The 
opacity broadening comes about because Mg II k is an op¬ 
tically thick line, with the peak separation of k 2 strongly 


Berger, T. E., De Pontieu, B., Schrijver, C. J., & Title, A. M. 
1999, ApJ, 519, L97 

Cauzzi, G., Reardon, K., Rutten, R. J., Tritschler, A., & 
Uitenbroek, H. 2009, A&A, 503, 577 
de la Cruz Rodriguez, J., De Pontieu, B., Carlsson, M., & 
Rouppe van der Voort, L. H. M. 2013, ApJ, 764, Lll 
De Pontieu, B., Berger, T. E., Schrijver, C. J., & Title, A. M. 
1999, Sol. Phys., 190, 419 

De Pontieu, B., Hansteen, V. H., Rouppe van der Voort, L., van 
Noort, M., & Carlsson, M. 2007, ApJ, 655, 624 


influenced by the chromospheric column mass. There is a 
significant sensitivity of the wing width to chromospheric 
temperatures (as well as microturbulence). 

The combination of these diagnostics with our numer¬ 
ical experiments indicates that the chromospheric tem¬ 
peratures in plage are remarkably constant and likely of 
the order of 6,000-6,500 K. This is compatible with the 
observed radiation temperatures of k 2 , and further con¬ 
firmed by the fact that the subordinate blend at 279.88 
nm is rarely found to be in emission, something that 
would be expected for higher temperatures. 

The nature of the significant microturbulence in plage 
remains unknown, but clearly the presence of strong 
shocks (driving dynamic fibrils and affecting line broad¬ 
ening; De Pontieu et al. 2015), torsional motions (De 
Pontieu et al. 2014a) in and around plage, and Alfven 
wave turbulence (van Ballegooijen et al. 2011) are can¬ 
didates. 

Currently ongoing numerical experiments suggest that 
the observed Mg Ilk profiles in plage can only be ex¬ 
plained by the combination of a TR at high column mass 
(i.e., hot, dense corona above) as well as strong chromo¬ 
spheric heating. Without the latter, densities at elevated 
chromospheric temperatures would be lower, thus leading 
to narrower profiles and a lack of filled-in Mg II k profiles. 

More work is required to better understand the cou¬ 
pling between chromospheric and coronal conditions in 
plage, as well as the apparent difficulty of reconciling the 
plage atmosphere considered here with observations in 
Ca II 854.2 nm and Ha, which may necessitate mod¬ 
els with multiple atmospheric components, perhaps sim¬ 
ilar to those found from modeling of Ca II H&K (Solanki 
et al. 1991) and Ca II 854.2 nm (de la Cruz Rodriguez 
et al. 2013). 
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